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Abstract: We study the finite- horizon optimal control problem with quadratic functionals for 
an established fluid-structure interaction model. The coupled PDE system under investigation 
comprises a parabolic (the fluid) and a hyperbolic (the solid) dynamics; the coupling occurs 
at the interface between the regions occupied by the fluid and the solid. We establish several 
trace regularity results for the fluid component of the system, which are then applied to show 
well-posedness of the Differential Riccati Equations arising in the optimization problem. This 
yields the feedback synthesis of the unique optimal control, under a very weak constraint on the 
observation operator; in particular, the present analysis allows general functionals, such as the 
integral of the natural energy of the physical system. Furthermore, this work confirms that the 
theory developed in Acquistapace et al. [Adv. Differential Equations, 2005] — crucially utilized 
here — encompasses widely differing PDE problems, from thermoelastic systems to models of 
acoustic-structure and, now, fiuid-structure interactions. 
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1 Introduction 

In this paper we consider the optimal control problem with quadratic functionals for 
a fluid-structure interaction model. Of major concern is well-posedness of the Riccati 
equations arising in the minimization problem, along with the feedback synthesis of 
the (unique) optimal control. The fluid-structure interaction is modeled by a system of 
coupled partial differential equations (PDE) comprising a Stokes system (the fluid) and 
a three-dimensional system of dynamic elasticity (the solid). The coupling occurs at 
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an interface separating two regions occupied, respectively, by the fluid and the sohd. It 
is worth mentioning at the outset that it is assumed that the motion of the sohd is due 
to inflnitesimal displacements. Accordingly, the fluid-solid interface is stationary; this 
and other modeling issues are discussed, e.g., in [17]. The mathematical description 
of the PDE system, that is the boundary control problem (j2.ip . as well as further 
literature will be given in the next Section. Our main goal is to establish the validity 
of a Riccati theory that would allow to control the structure, via boundary controls, 
acting as forces applied to the interface. 

It is well known that — even in the case of a single PDE — one of the main difficulties 
in a rigorous derivation of the feedback synthesis of the optimal control is the presence 
of boundary controls (or, more generally, unbounded control actions), combined with 
the lack of smoothing effects propagated by the dynamics (see, e.g., [10] and [23]). 
In fact, while the linear-quadratic control problem with unbounded control operator 
has a complete solution in the case of PDE models whose free dynamics is governed 
by an analytic semigroup, this solution may be out of reach in the case of other kind 
of dynamics. In particular, the case of purely hyperbolic PDE with boundary /point 
control is peculiarly differenlH from the parabolic case: it would suffice to recall that in 
finite time horizon problems the Riccati operator (or optimal cost operator) P{t) does 
not satisfy the differential Riccati equations, unless the observation operator possesses 
a suitable smoothing property. 

On the other hand, certain interconnected PDE systems combining parabolic and 
hyperbolic effects may give rise to an abstract control system y' = Ay + By which 
yields a singular estimate for the operator e^^B, near t = 0. This property — which 
is an intrinsic feature of control systems ruled by analytic semigroups — has been first 
identified in the analysis of an acoustic-structure interaction (where the overall semi- 
group was not analytic); see [3]. The essence of such estimates is the following: the 
parabolic component does induce a singular estimate (as a consequence of analyticity 
of the corresponding semigroup), while hyperbolicity 'transports' this estimate across 
the system through the coupling. Thus, if one can show that a singular estimate is valid 
for the entire system, then the theory in [2H [22] ensures a feedback control law with 
bounded (in the state space) gain operator, along with well-posed Riccati equations. 
This theory has been successfully applied to diverse composite PDE models, including 
some thermoelastic systems, beside to various acoustic-structure interactions. Several 
illustrations are contained in [22j and [24]; see also [15], |14j . and the recent |12j . 

For the fluid-structure interaction under investigation, which comprises a parabolic 
and a hyperbolic PDE, it was shown in [26] that a singular estimate (for the corre- 
sponding abstract evolution) is satisfled in the finite energy space, as long as the penal- 
ization in the quadratic functional does not involve the mechanical energy at a truly 
energy level. More precisely, the study in [26] established speciflc singular estimates 
and hence well-posedness of the Riccati equations in the special case of penalization 
of the mechanical variables below the energy level (say, sub-critial penalization), yet 

^ In the infinite time horizon case the so called gain (or feedback) operator B*P is intrinsically un- 
bounded and the analysis of the algebraic Riccati equations is subtle; see [T5], [33] . and the subsequent 
improvements in [32], |9]. 
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allowing full penalization of the fluid variable. 

The situation becomes much more difficult when the mechanical variables are pe- 
nalized at the critical level of the energy (see the functional (j2.9p l. In fact, not only 
the regularity results of [25j do not apply, but the theory pertaining to control systems 
which yield singular estimates ([M]) is no longer valid. (Indeed, if it were so, the 
gain operator would be bounded on the state space, while we will show that this is 
not the case; see Remark 12.81 ) 

The present work addresses the issue of solvability of the optimal control problem 
with general quadratic functionals (i.e. including critical penalization) for the PDE 
model ()2.ip . As we shall see, we provide (a positive) answer to the question remained 
open in [261, Remark 6.1]. This will follow in light of the theory introduced in [2], 
which is shown to cover the present case in view of the set of trace regularity results 
established and collected in Theorem 12.91 The theory developed in [2] is more effective 
in capturing the relevant properties of the dynamics, especially the ones which emanate 
from hyperbolicity. These ultimately allow to define the gain operator as an unbounded 
operator on a suitable functional space. In this respect, the variational aspect of 
the minimization process is critical in order to justify the arguments leading to well- 
posedness of Riccati equations. 

Let us recall that the optimal control theory in [2], while relaxing the 'singular 
estimate requirement', postulates other regularity conditions of global nature. This 
makes it possible to obtain meaningful solutions to the differential Riccati equations, 
despite the gain operator is not bounded on the state space. This, however, does 
not affect the synthesis, as the optimal solution still belongs to the domain of the 
gain operator. Originally arisen in the study of boundary control problems for an 
established system of thermoelasticity ([!]), so far this theory has been shown to apply 
as well in the case of certain acoustic-structure interaction model including thermal 
effects (Hg). 

The paper is organized as follows. In Section [2] we introduce the boundary control 
problem under investigation, along with the statements of our main results, namely 
Theorem 12.61 and Theorem 12.91 Moreover, we briefly record some necessary notation 
and the fundamental well-posedness result pertaining to the uncontrolled counterpart 
of the PDE system. Section [3] is entirely devoted to the proof of Theorem 12.91 which 
establishes the novel, distinct boundary regularity properties (of the solutions to the 
PDE system) which will ultimately result in solvability of the optimization problem, i.e. 
Theorem l2.61 Section S] contains the proof of Theorem l2.6^ based upon the application 
of the theory in [2]. Finally, a short Appendix collects the statements of the regularity 
results pertaining to the elastic component of the system — recently obtained in [7j and 
[26j — which are crucially utilized in the proof of Theorem 12. 9[ 

2 The PDE model, statement of main results 

The PDE model. The PDE model under investigation describes the interaction of 
a (very slow) viscous, incompressible fluid, with an elastic body in a three dimensional 
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bounded domain. Although the introduction of such models dates back to |27] . their 
PDE analysis has increased significantly only in the past decade. A mathematical 
description of the composite PDE system is given below. By and we denote the 
open smooth domains occupied by the fluid and the solid, respectively. Then C 
denotes the entire solid-fluid region, that is Q is the interior of U fi^. The boundary 
of is the interface between the fluid and the solid, and is denoted by Ts = d^ls- We 
finally denote by the outer boundary of 17/, namely = d^lj \dQs- It is assumed 
that the motion of the solid is entirely due to infinitesimal displacements, and hence 
that the interface Tg is fixed. 

The velocity field of the fluid is represented by a vector-valued function u, which 
satisfies a Stokes system in Qf, the scalar function p represents, as usual, the pressure. 
In the solid region the displacement w satisfies the equations of linear elasticity. 
(The density and the kinematic viscosity which usually appear in the Navier-Stokes 
equation are set equal to one, just to simplify the notation). The coupling takes place 
on the interface Tg. We recall from [T^ that the interface condition u = wt on Tg 
(in place of the usual no-slip boundary condition u = 0) accounts for the fact that 
although the displacement of the elastic body is small, its velocity is not (small, yet 
rapid oscillations). Thus, the PDE system is given by 



ut — dive(n) -|- S/p = 


in 


Qf 




X (o,r) 


div u = 
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Qf 
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In the above coupled PDE system, a and e denote the elastic stress tensor and the 
strain tensor, respectively, that are 

k=i ^* 

where A, fi are the Lame constants and 6ij is the Kronecker symbol. 

Since the present work is focused on the optimization problem, the subtle ques- 
tions related to the modeling of fluid-structure interaction phenomena, as well as to 
the analysis of well-posedness of the corresponding coupled PDE systems, will not be 
discussed here. Yet, well-posedness of the boundary value problem (j2.ip . with g = 
(that is the uncontrolled system (j2.3p below), is a prerequisite for the study of the 
associated optimal control problems. Thus, although many authors have contributed 
to the PDE analysis of nonlinear fluid-structure interaction models (where the dynam- 
ics of the fluid is ruled by a Navier-Stokes equation), existence of finite energy weak 
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solutions — even for the simpler Stokes-Lame system ()2.3p — has been an open question 
until O Theorem 2.2]. The reader is referred to [71 [8] for the analysis of well-posedness 
of the coupled PDE system (j2.3p : in addition, [8] includes a very clear introduction 
to the (nonlinear) fluid-structure interaction problem, along with a technical compar- 
ison with the previous mathematical literature. In Section 12.11 we shall recall, for the 
reader's convenience, the theory in [7] that is needed for our purposes. 

We finally note that while the present study follows the variational approach of 
[7], exploiting the novel boundary regularity results established therein, semigroup 
well-posedness and stability properties of the linear model have been investigated in 
[1]; see the survey paper [5j and its references. For the uniform stabilization problem, 
see [6]. 

Further references. There is a large literature on coupled fluid-structure evolution 
problems. Most works address the issue of developing models for specific physical 
problems and/or their numerical simulation. Two main different scenarios arise from 
the applications: the case in which the fluid is flowing in a tube with elastic walls, 
such as the blood through arteries, and the case where one or more elastic bodies are 
immersed in a fluid flow. The PDE model under investigation pertains to a physical 
situation falling under the latter category. 

A very nice introduction to fluid-structure interaction problems is provided by [T7] . 
Recent treatises with focus on modeling and numerical analysis are |30j and [29]. An 
in-depth PDE analysis of well-posedness of these nonlinear models has indeed appeared 
only recently. Relevant contributions to this problem are given (without any claim of 
completeness) by [M], [H], the aforesaid [IT], [l6], [11], [3 [8] and, lastly, [20]. For 
more information on this subject, see the bibliography therein. 



2.1 Variational and semigroup formulation 

Before giving the statement of our main results, let us preliminarly recall from [7] 
some basic notation, and the chief facts which pertain to the uncontrolled problem, 
that is system (12. ID with g = 0. Further technical results obtained in [7J and [26j will 
be needed in the proof of our main result; these will be recorded in an Appendix for 
convenience. 

The uncontrolled model. Let us introduce the free system corresponding to (j2.ip . 
namely 
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(2.3) 
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The energy space for the PDE problem (j2.3p is 

Y = H X H^Qs) >^ L2ins) , 

where H is defined as follows: 

H := {u e L2{^lf) : divn = , u • z^|r^ = O} . 

In addition, we denote by V the space defined as follows: 

V := {v e H^{nf) : divu = , u|r^, = O} ; 

we shall use the following distinct notation for the various inner products which will 
occurr throughout the paper: 



{u,v)f := / uvd^lf, {u,v)s ■= / uvdil.s, {u,v) := / uvdTs- 
The space V is topologized with respect to the inner product given by 

{u,v)ij := / e{u)e{v)dClf ; 



the corresponding (induced) norm | • \ij is equivalent to the usual ff^(r2j) norm, in 
view of Korn inequality and the Poincare inequality. 

Remark 2.1. The norm || • in the Sobolev space H'^{D) will be shortly denoted 

by I • \r,D throughout the paper. Note that all the Sobolev spaces related to u and 
w are actually (H''')^: the exponent is omitted just for the sake of simplicity. 

Let us recall from [7] the definition of weak solutions to the (uncontrolled) PDE 
system ()2.3p . 



Definition 2.2 (Weak solution). Let (uq, wqjWi) € H and T > 0. We say that a triple 
(n, w, Wt) € C([0, T],H X H^^Qg) x L2{^s)) is a weak solution to the PDE system (|2.ip 

• {u{-,0),w{-,0),wt{-,0)) = {uo,wo,wi), 
. ueL2{0,T;V), 

. a{w) ■ V G L2(0,r;//-V2(r,))^ = u\v^ E L2(0, T; /7i/2(r^))^ 

• the following variational system holds a.e. in t ^ {0,T): 

(i{u,(l))f + (e(u), e{cj)))f - {a{w) .„ + g,cP) = 



(2.4) 

j^{wt,i))s + {(T{w),e{'4>))s - {(t{w) • = 



for all test functions (p & V and ip € H^{i^s 
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Remark 2.3. It is important to emphasize that the regularity properties of the normal 
stresses (see the third item of Definition 12. 2p do not follow from the interior regularity 
of the fluid-structure variables. It is an independent regularity result, showing the ex- 
ceptional behavior of hyperbolic traces. This regularity property is necessary in order 
to justify the variational definition of weak solutions (see (j2.4p ). While there are other 
definitions of solutions to nonlinear PDE models of fluid-structure interactions which 
do not require additional regularity on the boundary (see, e.g., [27], [IT], [1]), yet these 
definitions are not adequate to variationally decouple the (finite energy) weak solu- 
tions of the two equations. On the other hand, this decoupling is crucially important 
in the present analysis, aimed at identifying the distinctive regularity properties of the 
overall dynamics, that play a major role in the study of the associated optimal con- 
trol problems. Exploiting the distinct features (analyticity and hyperbolicity) of the 
decoupled dynamics makes it possible to establish the sharpest results for the coupled 
PDE system. (This fact was recently utilized in [20], as well.) Consequently, the issue 
of "hidden" regularity of the hyperbolic component is central to the problem studied 
and its solution. 

Existence of weak (global) solutions of a nonlinear generalization of the PDE prob- 
lem ()2.3p has been established in [7j. 

Theorem 2.4 (Existence of weak solutions, [7J). Given any initial datum (uq, wo,wi) £ 
Y and any T > 0, there exists a weak solution {u,w,wt) to the system ()2.3p such that 

VwL GL2(0,T;/7-V2(r,)), = WtL G L2(0, T; //^/^(r,)) . 

The control system, semigroup formulation. Aiming to apply the optimal 
control theory pertaining to a general class of evolutions — in the present case, the one 
developed in [2] — it is convenient to recast the boundary value problem (|2.ip as an 
abstract control system in a Hilbert space. Accordingly, let us introduce the fluid 
dynamic operator A : V ^ V', defined by 

iAu,(t>) = -{e{u),e{4>)) y^eV, 

and the (Neumann) map N : L2(Ts) — > H defined as follows: 

Ng = h^{e{h),e{<f>)) = {g,^) y<peV. 

The chief properties of the operators A and are recalled in the Appendix. 

Then, if we set y = {u,w,wt), the boundary value problem (j2.ip reduces to the 
linear control system 

(y' = Ay + Bg in [V{A*)]' 5) 

\y{o) = yo 

where the (dynamic) operator A : T^{A) C y ^ y is defined by 

A ANa{ ■)-u 
0/ 
divo-(-) 



A AN(7{ ■)-v 
^ = I / ) , (2.6) 
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with domain 



V{A) = {y = {u,w, z) e H : u e V, A{u + Na{w) ■ u) e H , z eH^i^s 
divaiw) G L2(0s) , z\r, = u\r,} , 

and the (control) operator B : L2{Ts) = U [V{A)]' reads as 




(2.7) 



Given a quadratic functional 



Jig)= r {\ny{t)\l + \g{t)\l)dt, (2.8) 
Jo 

the optimal linear-quadratic (LQ) control problem is to minimize the functional (|2.8p . 
over all control functions g E L2(0, T; [/), with y solution to (|2.5|) corresponding to g. 
As already pointed out in the Introduction, we aim to include in the present analysis 
non-smoothing observation operators TZ, such as the identity operator; hence, TZ is 
initially assumed to satisfy just TZ G C{Y, Z). By doing so we admit natural quadratic 
functionals such as the following. 

Jig) = \l { \^it)\o,n, + Kw{t)), eiwms + \wtit)\ln, + \g{t)\lr, } dt (2.9) 
which penalizes the full quadratic energy E{t) of the system. 



Remark 2.5. We already emphasized that the study performed in [26] did not pro- 
vide solvability of optimal control problems with general quadratic functionals: in 
particular, it did not cover the case of natural functionals such as ()2.9p . On the other 
hand, the analysis carried out in [2^ — despite the final constraint on the observation 
operator TZ — included the case of Bolza problems, where the penalization affects also 
the state at the final time T < oo, namely, when the functional to be minimized is 
given by 

J{g)= [ {\TZy{t)\l + \g{t)\l)dt + {gy{T),y{T))w. (2.10) 
Jo 

Note that the LQ-problem with Bolza-type quadratical functionals is not discussed 
here. In fact, the LQ-problem with quadratic functionals of the form (j2.10p (with 
Q ^ 0) for the class of control systems (12. 5j) described by the Assumptions 14.11 has 
not been investigated yet. 



2.2 Statement of the main results 

The main result of the present work is the proof of well-posedness of the (differential) 
Riccati equations corresponding to the optimal control problems associated with the 
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fluid-structure model (j2.ip . along with all the inherent assertions about solvability of 
the optimization problem; see Theorem 12. 6i This variational result, however, critically 
relies on the novel trace regularity results established specifically for the (uncontrolled) 
PDE system (|2.3p in Theorem 12.91 which thus constitute the major technical contri- 
bution of the present work. As we shall see, the proof of this set of regularity results 
is based on the interplay between the maximal parabolic regularity of the fluid com- 
ponent with the 'hidden' regularity of the traces of the hyperbolic (solid) component. 
Indeed, the fact that the coupling is of hyperbolic/parabolic type will be critically 
utilized. 

2.2.1 The solution to the optimization problem 

With reference to the PDE model introduced in the previous section, let us consider 
the optimal control problem ()2.5p ~ (|2.8p . that is 

Minimize the functional J{g) in (j2.8p . over all g G ^2(0, T; L2{Ts)), where y(-) = 
y{''-,VQ-,9) solves the control system (|2.5p . 

Then we have the following. 

Theorem 2.6. Consider the optimal control problem (j2.5p ~ (j2.8p . with A and B given 
by (|2.6p and (j2.7p . respectively. If the observation operator satisfies 

n*n G £{V{A'),V{A*')) (2.11) 

for some e G (0, 1/4), then the following assertions hold true. 

1. For any initial state yQ gY there exists a unique optimal control g^{-) G -^2(0, T; L2{Ts 
such that 

J(g°) = min J(g) . 

gGL2(0,T;L2(r.) 

The optimal pair {g^ {■) , {■)) has the following additional regularity: 
yO(-) = [u^i-),w^{-),w',{-)] G C{[0,T];Hx HHn,) x 12^); 
5°(-)G n Lp{0,T;L2{Ts)). (2.12) 

l<p<oo 

2. There exists a non-negative, self adjoint operator (the Riccati operator) P{t) G 
^{Y), t G [0,T], defined explicitly in terms of the data, such that 

J(5°) = (P(0)yo,yo)y; 
more precisely, P(-) G C{Y,C{[0,T],Y)). 

3. The gain operator B*P{-) satisfies B*P{-) G C{V{A^),C{[Q,T],L2{Ts))); more- 
over, one has (the feedback synthesis of the optimal control) a.e. in [0,T]: 

g\t) = -B*P{t)y\t), VyoGl^. (2.13) 
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4- The operator P{t) is the unique solution of the Differential Riccati Equation 
satisfied for < t < T and x = {xi,X2,x^),y = (yi, 1/2, ys) ^ D{A), 

j^{P{t)x,y)Y+{-A*P{t)x,y)Y+{P{t)Ax,y)Y+{T^x,ny)Y = {B*P{t)x,B*P{t)y) , 

(2.14) 

with 

lim (P(t)x,x) =0 Vx G y . 

Remark 2.7. Since here the dynamics operator A is the generator of a s.c. contraction 
semigroup with A~^ € C{Y), then the domains of fractional powers 'D{A'') in (j2.1ip 
may be computed as intermediate spaces between 'D(A) and Y. The same holds for 
V^A*"^). (For a comprehensive hst of cases where the identity ['D{A),Y]i-0 = 'D{A^) 
holds true, see, e.g., [23^ § 0.2.1]). Then, it is not difficult to show that in the present 
case 'D{A*'^) = T>{A'^), provided e is sufficiently small. Therefore, assumption (j2.1ip 
is satisfied, with a non-smoothing observation operator, such as TZ = I. This natural 
situation was indeed left as an open problem in |26j . 

Remark 2.8. Observe that the optimal pair does not display the typical regularity (in 
time) exhibited in the case of control systems whose underlying semigroup is analytic 
(or, more in general, when singular estimates are satisfied). In particular, the optimal 
control is not continuous. This is not surprising, in view of the influence of both 
hyperbolic and parabolic effects on the overall behavior of the solutions. 

Moreover, the gain operator B*P{t) is no longer bounded on the state space Y , 
but just densely defined. However, this does not affect the final result, as the feedback 
formula holds for any initial state in the finite energy space. Thus, the observation 
operator 7^ need not have regularizing effects, and IZ can be critical. 

The above observations also explain why the previous Riccati theories are intrinsin- 
cally unapplicable in the critical case, as they lead to hounded gain operators, in con- 
trast with the case under examination. 



2.2.2 Trace estimates 

The fundamental analytic tool which will enable us to show Theorem 12.61 is a complex 
of boundary regularity results pertaining to the fuid component of the PDE problem 
(|2.3p . These traces' regularity estimates of u (and ut) on the interface are the 
"PDE counterpart" of the abstract regularity properties of the (unbounded) operator 
^*g.A*t needed to invoke the optimal control theory of [2]. These regularity estimates 
are, however, also of independent interest. 

Theorem 2.9 (Traces' regularity). Consider the uncontrolled Stokes-Lame system, 
namely the PDE system ()2.3|) . Let y{t) = {u{t),w{t),wt{t)) be the solution corre- 
sponding to initial data yo = {uq,wq,wi). Then the fluid component u admits a de- 
composition u{t) = ui{t) + U2{t), and the following statements pertain to the regularity 
of the traces of ui, U2 and ut on Tg, respectively. 
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(i) The component ui satisfies a pointwise (in time) "singular estimate", namely 
there exists a positive constant Ct such that 

Mt)\L,[T.)<-^MY yyoeY, vt€(o,r] (2.15) 

(with arbitrarily small 5 > 0). 
(a) The component U2 satisfies the following regularity: 

(iia) if yo G Y, then U2\r^ G Lp{0,T; L2{Ts)) for all (finite) p > 1; 

(iib) tfyo G V{A'), e G (0, i), then U2\r., G C([0, T], L2(r,)). 

(Hi) Let now yo G V{A^^^), with 9 G (0, |). Then, the fluid component u of corre- 
sponding solution satisfies, for some q G (1,2), 

ut\r^ eLg{0,T-L2iTs)) (2.16) 
continuously with respect to yo, that is there exists a constant Ct such that 

\\ut\\Lg{0,T;L2{rs)) ^ C't ||yO • (2-17) 

The exponent q will depend on 0: more precisely, given 6 G (0, |), one has 

1<.<J^. (2.18) 

The remainder of the paper is devoted to the proof of the two main results stated 
in Theorem 12.61 and Theorem 12.91 Section [3] deals with the above boundary regularity 
results, which will be next utilized in Section H] to establish Theorem | 



3 Proof of the trace regularity results 

This section is entirely devoted to the proof of our main contribution, that is Theo- 
rem E^l 

Proof of Theorem \2.9\ Our starting point is the equation satisfied by u{-), namely 
ut = Au + ANa{w) ■ whose solutions are given by 

u{t) = e^*uo + f e^^^-'^ANa{w){s, ■)u ds ; (3.1) 
Jo 

the above expression yields the natural splitting u{t) = ui{t) + U2{t), with 

ui(t) := e^*uo, U2{t):=[ e^^^''^ ANa{w){s, ■)u ds . (3.2) 

Jo 
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In view of N*Au = —u\r^ (see Lemma lA.ip . the corresponding traces on Tg read as 
u^\^^ = -N*Aui(t) = -N*Ae^^uo , 

nsir, = -N*Au2{t) = -N*A I e"^^^''^ ANa{w){s, -yds , (3.3) 



respectively. 

( i ) The singular estimate in (j2.15p follows as an immediate consequence of the well 
known estimates pertaining to analytic semigroups: 



uo\ 



(3.4) 

< \\N*A^I^-'\\ \A^'^+'e^'uo\Y ~ CTt-^''~'\uo\ . 
This shows the validity of assertion (i). 

(ii) Let initially yo = {uq,wq,wi) & Y. In view of (|3.3p . it is clear that the sharp 
regularity theory pertaining to the the wave component will play a central role in the 
study of the regularity of the traces of U2{t) on Tg. More precisely, we shall utilize 
the recent trace results obtained in [71 Theorem 3.3] and refined in [26', Lemma 5.2]; 
see Lemma lA.2l in the Appendix. Accordingly, following the decomposition of a{w) ■ v 
established in Lemma lA.2( it is convenient to introduce a further splitting, namely 

U2(t) = / e^(*-^UiVa(u'i)(s,>(is+ / e^^^-'^ AN a{w2){,s,-)v ds . (3.5) 







1*21 (t) 1*22 (t) 

Thus, one has first 



N*Au2x{t) = N*A / e^^^~'^ANa{wi){s) ■ v ds 
Jo 

= [Ar*A3/4-]Ai/4+^+i/2 re^(*-^) A'/^Na{wi){s) ■ v ds 





fis) 

where / € C{[0,T], L2{0,s)) in view of Lemma I A. 2 1 and Lemma lA.ll Consequently, 

[\^^'-'^f{s)ds € C{[0,T],D{A'-n), (3.6) 
Jo 

with arbitrarily small cr > 0; see, e.g., [231 Proposition 0.1, p. 4]. Therefore N*Au2i G 
C([0, T], L2{Ts)), and a fortiori we obtain 

N*Au2i G Lp(0, T;U) Vp > 1 . (3.7) 

As for the second summand N*Au22it), because of the different regularity of 
a{w2){s) ■ V we rewrite in a different fashion: 

A^Mu22(t) = \N*A^I^-'\ A^l^^'^' /"*e^(*-^) {A^I^-'Na{w2){s) ■ v\ ds . 

Jo ' V ' 
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Notice now that the above integral is the convolution 

rt 

K{t — s)^p[s) ds , 



with € L2{0,T;U) and the kernel K such that ||ir(s)|| ~ ^7^+271 where e can be 
taken arbitrarily small. Hence K G L2_o-(0, T; [/) for arbitrarily small a > 0. Thus, 
the Young inequality yields 

N*Au22 G Lp{0, T-U) Vp > 1 . (3.8) 

Thus, ()3.8p combined with (13. 7p shows the validity of assertion (na). 
Let now uq € D{A'')^ e > 0. In this case we have 

and by Lemma [A. 3 1 it follows u\r^ G H'^{T,s), provided that e < |. This enables us to 
apply the second part of Lemma IA.21 which gives 

a{wi) ■ V G C([0, T], i/-i/2(r,)) , a{w2) ■ v e . (3.9) 

Now, the analysis of N* Au2\ follows closely the one in item (iia), yielding the conclu- 
sion Au2\ € C([0, r], L2(rs)) (this is justified by the membership (j3.6p ). Instead, 
on the basis of the novel regularity of 0(^2) • v in (|3.9p . from parabolic theory it follows 
that 

SO that 

N''Au22 G F^+i'^/2+i/2(5.^) ^ FV2+^/2(o,r;L2(r.)) C C([0,T],L2(r,)). 

As both iVMn2i and N* AU22 belong to C([0, T], L2(r,)), then N* Au2 G C([0, T], L2(r,)) 
and {iih) is proved. 

(Hi) In this last step we aim to ascertain the regularity of the boundary traces of the 
time derivative ut on Y g. We return to the mild solution (j3.ip and compute 



ut{t) = Ae^^ + A e^^^-'^ANa{w){s, ■)v ds + ANa{w){t, >^ (3.10) 



fi(i) 

V2it) 

which can also be rewritten as 

V2l{t) 



V22{t) 



ut{t) = Ae^\iQ + j e'^^*~''^ANcj{ws){s,-)yds + Ae^^Na{w){<d,-)u. (3.11) 

Vl{t) ^ V ' 

V2(t) 
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The plan we aim to carry out is to discuss first the regularity of the function V2 ■= -gf 
when yo S D{A) by using its expression in (j3.1ip . Next, when ^ Y , we would 
rather utilize (j3.10p . and then use interpolation arguments to establish the regularity 
corresponding to initial data in D{A^~''). Only subsequently we shall derive the trace 
regularity of V2 by applying the operator —N*A. 

When yo ^ D{A), by standard semigroup arguments we known that cr{wt) • v 
exhibits the same regularity as that of a{w) ■ v when yo ^ Y , i.e. (invoking once again 
Lemma IA.2P 

a{wt) • = ai + ^2 G C([0, T], ii-i/2(r,)) © ^2(0, T; L2(r,)) . 

To pinpoint the regularity of V21, we now utilize the above splitting and follow the 
analysis carried out in the proof of {ii). More precisely, combining elliptic regularity 
(in particular, Lemma lA.ip . with the analyticity of the semigroup e"^*, along with the 
(singular) estimates pertaining to A^e"^^, we first obtain, for any t and any 5 < 1/2, 

1^1/2-5 /"*e^(*-^)^iVc7i(s)ds| = \A^~^ f e^^^-''>A^/^Nai{s)ds\ 
Jo Jo 

/■* 1 

(t - sy-^^^ ll^illc([o,T],//-i/2(r,)) < C'lkillc([o,T],/f-i/2(r,)) • (3-12) 



As for the latter term, we apply as well A^^^ ^ and rewrite as follows: 

^1/2-5 I" ^ f\^3/4-5/2^Ait-s)^ [^3/4-^2^] ^3(5) ds , (3.13) 

Jo Jo 

where it is clear now that the integral is the convolution of and L2 (in time) 
functions, respectively. On the strength of the Young's inequality, we get L4-regularity 
in time, so that 

V2i G C([0, T],D{A'/^-^)) e L4(0, T; D{A'/^-^)) , < <^ < ^ . 

This implies the membership 

V2i€L^{0,T;D{A'^^-^)), 0<5<^. (3.14) 

On the other hand, still with yo G D{A), one has just a{w) ■ v G C([0, T], iJ-i/2(rs)) 
which suggests us to rewrite ^22 as follows: 

V22{t) = A^l^e^' (^i/2iVa(«;)(0, » ; 

then, again in view of Lemma lA.ll and of the usual singular estimates pertaining to 
analytic semigroups, it follows 

^22 GLg(0,r;D(^i/2-^)), (3.15) 
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provided that q{l — 6) < 1. Therefore, ()3.14p combined with ()3.15p yields, for any 
< (5i < i 

2/0 G D{A) ^V2£ L,,(0,r;Z)(Ai/2-5i)) = Lg,{0,T; H^-^^' {nf)) , (3.16) 
where qi E (1,2) depends on 5i; more precisely, 

1- di 

Let now yo G Y. In this case we use the decomposition (j3.10p . and begin with the 
analysis of V2a- Setting w = wi + W2 (according with Lemma [A. 2p . one has 

A f e'^^^-'^ANa{wi){s,-)vds = 
Jo 

f e^^'-'^A^/^Na{wi){s,-)iyds e C([0,T], [D{A^/^+'^)]') , 
Jo 

while 

a[ e^^*-''>ANa{w2){s,-)iyds = 
Jo 

A^^+'^A f e'^'^'''^A^/^-'^Na{w2){s,->ds e L2{0,T;[D{A^/^+'^)]') 
Jo 

where both ei and €2 can be taken arbitrarily small. As a result, 

V2a G L2{0, T; [D{A'/^+^)]') , < e < ^ . (3.18) 
As for the term V2b, readily 

ANaiwi){t, ■)u = A^/^A^/^Na{wi)it, ■)v G C([0, T], {D{A^I'^)\) 

while 

ANa(w2){t, ■> = A^'^^^[A'l^-^N]a{w2){t, ■)v G L2(0,r; \D{A^l^^^y{) , 
and since e can be taken arbitrarily small, we deduce as well 

V2b = ANa{w){t, ■)u G L2(0, T; [D{A'/^)]') . (3.19) 
On the basis of ()3.18p and ()3.19p . we obtain 

yoeY^V2 = V2a + V2b S i>2(0,r; [Z)(Ai/2+fc)]') = L2(0,r; (17^)]') , (3.20) 

if < (52 < i. 

Thus, (j3.20p . combined with (j3.16p . gives by interpolation 

yo(^D{A^~^)^V2eLg,{{),T;W), (3.21) 
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where qi is as in ()3.17p and W is the interpolation space 
if 

s = (1 - 2(5i)(l -9)- e{l + 262) = 1 - 2(5i - 20(1 + (52 - 5i) > ; 

see {28^ Theorem 12.5]. Notice that by taking, for instance, 5i = 82 ='■ S, one has 
s > 1/2 provided that 

e + 6<^. (3.22) 

In this case V2 G H^{0,f) with s > 1/2 and hence its trace on Tg is well defined. Notice 
that, in view of the constraint (j3.22p . we need to require < < j. Consequently, 
given any such that < < ^, choosing, e.g., 5 = 1/4 — in view of (j3.22p . from 
(fOT]) it follows 

yo G D{A'-')) =^ N*Av2 G Lg,{0,T; L2{T,)) Vgi < . (3.23) 

It remains to establish the regularity of the first summand N*Avi{t) = N*Ae^^Auo 
when yo G D{A^-'^). In this case Uo G (F^(0/), L2(0/))e = H^-^Vtf), and from 

N*Avi{t) := N*Ae^'Auo = [N* A"^-']A^'^+'+^'^+'I\^' A^^-'^'^uq , 
it immediately follows 

yo G D{A^-'') =^ N*Avi G L,,(0, T; L2(r,)) Vg2 < j^^^0 + Ae ' ^^'^^^ 

Notice that in the above membership the Sobolev exponent q2 belongs to (1,2), as 
well. In conclusion, since e in p.24p can be taken arbitrarily small, the regularity 
Lq^{Q,T;L2{Ts)) combined with Lq^{Q,T-L2{Ts)) (in (l3:23D and IJ^TM . respectively) 
imply the membership 

yoeD{A^-')^ut\T^=:v\r^^Lq{Q,T-L2{Vs)) < ^-^^ , (3.25) 

which concludes the proof. 

4 Proof of Theorem 12.61 

The conclusions stated in Theorem 12 . 6 1 will follow from [TJ Theorem 2.3], once we verify 
the standing assumptions, which are recorded below for the reader's convenience. 

Assumptions 4.1. For each t G [0, T], the operator B*e^^ can be represented as 

B*e-^*'yo = F{t)yo + G{t)yo, t>0, yo e V{A*) , (4.1) 

where F{t) : Y ^ U and G{t) : ^{A*) ^ U , t > 0, are bounded linear operators 
satisfying the following assumptions: 
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1. there is e (^,1) such that \\F{t)\\c(Y,U) < C t"^ for all t G (0,T]; 

2. the operator G{-) belongs to £.{Y, L'''(0,T;U)) for all p € [l,oo), with 

\\G{-)\\c{Y,Lp{o,T;U)) < Cp < CO Vp e [l,oo) ; (4.2) 

3. there is e > such that: 

(a) the operator G{-)A*''' belongs to C([0, T], [/)), with 

\\^~'G{ty\\c(u,Y) < c < oo Vt G [0, T] ; (4.3) 

(b) the operator 71*11 belongs to C{V{A'^),V{A*^)); 

(c) there is q G (1,2) (depending, in general, on e) such that the operator 
B*e-^' n*nA^ has an extension belonging to C{Y, Li{0,T;U)). 

Remark 4.2. Note that the set of requirements in Assumptions 14.11 involves the 
regularity (in time) of the operator jB*e"^**, both locally at the origin (with singularity 
controlled by 7), and globally (in Lp sense). 

We now prove that the regularity results established in Theorem 12.91 implv all the 
requirements in Assumptions 14.11 with suitable values of 7, e and q. 

Verification of Assumptions \4-l\ Following |26t Proof of Theorem 5.1], it is not 
difficult to verify that given any initial state yo = {uq,wo,wi) G 2?(^*), one has 
B*e-^**yQ = N*Au{t) = —u{t)\r^, where u{t) is the first component of the solution 
y := {u,w,wt) to the (homogeneous) adjoint system 

y'{t) = A*m 

y(o) = yo ■ 

The abstract adjoint system yields a system of coupled PDE which is essentially the 
same as the original boundary value problem ()2.3p . except for few changes of sign in 
the equations. As a consequence, the regularity results established by Lemma IA.2I 
and Lemma [A. 3 1 readily produce analogues, by replacing 'D{A) and y = {u,w,wt) by 
'D{A*) and y = {u, w, wt), respectively. Similarly, the fluid component u of the solution 
y to the dual PDE system satisfies — mutatis mutandis — the regularity properties in 
Theorem 12.91 

1. In light of the decomposition of n found in Theorem 12.91 let us set 

F{t)yo := ui{t)\ , G{t)yo := U2(t)L . 

Then, the first statement in Theorem 12.91 along with the estimate (j2.15p . provides us 
with the soughtafter singular estimate, with (optimal) exponent 7 = 1/4 + 5, and the 
first of Assumptions St] is satisfied. 

2. Assertion (iia) in Theorem 12.91 is nothing but the regularity condition 2. of the 
Assumptions HTH valid for all p G [1, 00). 
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3. Condition (iib) of Theorem [2J] translates into G{t)A*~'yo G C([0, T], L2(r,)), 
which in turn gives the assertion (3a) of the Assumptions 14.11 with no constraints 
on e. It remains to verify the tricky assertion (3c) of Assumptions 14.11 This will be 
implied by condition (in) in Theorem 12 .91 We first claim that the estimate (j2.17p in 
(Hi) of Theorem 12.91 yields, for any 6 G (0, 1/4), the following one: 

\\B*e^'U*%o\\L'>io,T;L,(rs))<C\yo\Y, yoeV{A*'), i < q < , (4.4) 

This is easily seen if one just observes that if yo G ^{A^^) one has 

B e-^ A yo = B e-^ A [A yo) = B —e-^ zo = tit|r, , 

where now zo := A*^'^yo G V{A*^-^). Then, in view of the assumption (3b) on the 
observation operator TZ, one concludes 

||0*e■^**7^7^M^^/o||L,(o,T;L.(^.)) = ll^*e-4*U*'X-'7^7^M^o||Lno,T;L.(^.)) 

<C\yo\Y, yoG^(^'), 

i.e. condition (3c) is satisfied with e = 0, < < 1/4, for any q such that 1 < q < 
4/(3 + 4^). This completes the proof of Theorem \2M □ 

A Appendix 

For completeness' sake and for the reader's convenience we record some results shown 
in [7] and in [26], which are used frequently or critically in the proof of our main result. 

Lemma A.l (|26j). The Green (Neumann) map N : L2{Ts) H satisfies the follow- 
ing regularity results. 

(i) One has N*Au = — nlr^; u gV, where the adjoint is computed with respect to the 
L2-topology. 

{ii) N e £(L2(r,),P(A3/4-5)) n/:(i7-i/2(r,),P(^i/2)) for any5,^<5< f . 

Lemma A. 2 (0123). Let {wq,wi) E iJ"+^(0^) x i?"(0^), with < a < \, and let 
f G L2(0, T; //^/^(r^)). Then, the solution of the initial/boundary value problem 

Wtt — div(T(w) =0 in Qs 

< j-^w\v, =f onlls (A.l) 

^'u;(0, •) = ■Wo , wt{^, ■) = wi in 9.s 

can be decomposed as w = wi + W2, where a{wi) ■ v G C([0, T], //^^/^(T^)), while 
a{w2) • y G L2(0, T; L2(rs)). If, in addition, f G H°'{'Es), then a{w2) • v G -ff"(Ss). 
Moreover, the following estimates hold true. 

hiwi) ■ '^llc([o,T],H-i/2(r,)) < C'i(iu^o|?,n, + ki|o,n, + l/lL2(o,T;//V2(r,))) 
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Lemma A. 3 ([M])' Consider the uncontrolled counterpart of the PDE problem (j2.ip . 
that is (j2.ip with g = 0. Let initial data satisfy {uq,wo,wi) € L2{0,f) x x 
< a < Then, for any T < oo we have u\r^ € H°'{'Es) and the following 
estimate holds true: 

Remark A. 4. It it important to emphasize that the original proof of Lemma IA.3I 
(given in \26\ Lemma 5.3]) estabhshed local (in time) regularity for any T < Tq, with 
sufficiently small Tq. However, a slight variant of the proof yields the estimate (|A.2p 
without any constraint on T (provided it is finite). Indeed, it is sufficient to observe 
that the inequality (58) in the original proof, that is 

|t^3k"(s,) < K (|n|^^a(2,)T^+2^ + +\wo\i+a,ns + kiU.f^J , (A.3) 
can be actually improved to 

|t^3|//-(s,) < ekU"(E,) + C'el^iUaCs.) + {\wo\i+a,ns + \wi\a,ns) ■ 
The above estimate follows as a consequence of the following interpolation inequality, 

£ \D?-'u\l(^rJt < e£ |A"^lL(r.)rfi + CM^r.^dt , 

which makes it possible to avoid the use of T^/p in order to control the size of the 
principal term in the estimates preceding ()A.3p . 
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